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Abstract 

Formulas are derived to compute the mean number of times a site has been visited during 
symmetric Levy flights. Unrestricted Levy flights are considered first, for lattices of any dimension: 
conditions for the existence of finite asymptotic maps of the visits over the lattice are analysed 
and a connection is made with the transience of the flight. In particular it is shown that flights 
on lattices of dimension greater than one are always transient. For an interval with absorbing 
boundaries the mean number of visits reaches stationary values, which are computed by means of 
numerical and analytical methods; comparisons with Monte Carlo simulations are also presented. 
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I. INTRODUCTION 



Levy flights are a model of diffusion in which the probability of a |z|-length jump is 
"broad", in that, asymptotically, p(z) ~ |2;|~ 1_a , < a < 2. In this case the sum Xk = ^ Z{ 
is distributed according to a Levy distribution, whereas for a > 2 normal diffusion takes 
place 

Q,0 

. Interesting problems arise in the theory of Levy flights when considering the 
statistics of the visits to the sites, such for instance the number of different sites visited 
during a flight 

0,0; 

in this paper we consider a different, but related, problem, namely 
the number of times a site visited by a random flyer. 

Suppose that a random walk takes place on a <i-dimensional lattice £, let s be a site of 
£ and let Pu (s) be the probability that after k steps the walker is at s. The mean value of 
visits to the site s after n steps is p 

n 

Mf(s)=5>f(s); (1) 

fc=0 

n 

since derivation of Eq. (0) does not depend on the specific form of the walk p , it holds 
also for Levy flights. In the following it will be assumed Mq(s) = Pq(s) = <5 Sj o; the asvmptotic 
value of Mn \ denoted by M^ d \ is defined as M.^ = lim^oo M^f 1 . It is known [6( that a 
random walk is transient if and only if ( s ) < 00 ! ^ n other words the existence of 

finite Ad^ implies that the walk is transient. 

Levy flights have a wide range of applications (see for instance |7| and references therein) 
and, in particular, analysis of the number of times a site is visited can be relevant in those 
processes, such as random searches, in which it is important not just to determine what sites 
have been visited but how often they have been visited; examples of possible applications 
range from animal foraging to exploration of visual space 0. Moreover can be 

given the following interpretation, useful for possible applications: assume that particles 
undergoing a Levy flight are continuously generated at the origin, then, at time n, (s) oc 
C { n\s), where c£ d) (s) is the mean number of particles at site s ^ |l0|. This property of 
Mn has been used, in a model based on electrons Brownian motion, to simulate distributions 
ofe missi v ity „f S „_0. 
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II. INFINITE LATTICES 



Consider first one-dimensional, infinite lattices; the probability of occupancy of site x 



after k steps is 



111 



P$ii*)= E v{x-y)p£\vl (2) 



where p(y) is the probability of having a displacement of y sites. In case of symmetric 
Levy flights the canonical representation of p and are 



,w g Q 



1 f°° 

P{y) = - cos qyexp(-cq a )dq, (3) 

i r°° 

Pu^(x) = — cosgxexp(— ckq a )dq, (4) 
7T Jo 

where < a < 2 and c is a real number, which in the following will be set equal to 1 for 
simplicity Q]; a scaling relation holds between pjp and o, namely pf?\x) = k^ l / OL p{xk^ 1 ^ a ). 
If a = 2 Eqs. (j3J), (jlj) yield the Gaussian distribution Lj], jil, whereas, if a > 2, P^ 1 fails to 
be a proper distribution not concentrated at a point 6]; therefore representations (j3J), 
are valid only in the interval < a < 2. More recently it has been shown that the analytic 



forms of p and P^ can be given through a Fox function |12l 

Application of ((H) and of the scaling relation leads to Mn\x) = 8 Xi0 + YHl=i k~^p {-jjjz), 
and in particular, recalling that p(0) = (7ra:) _1 T(l/a;) Q], 

M«(0) = l + ^5>-i, (5) 

fc=i 

with X)fc=i converging to a finite value for n — ► oo if and only if a < 1 [l3j; in this 

Al«(0) = l + ^Mc(l/«), (6) 



case 



7T« 



where £ is the well known Riemann zeta function [13]. Thus Eqs. © and (JBJ), show that 
the visit to site x = is a transient state if and only if a < 1. 

The trend of M^\o) as n increases can be computed by making use of the formulas 
related to the zeta function Q|; for a < 1 the result is 
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M W (0 ) = l + EeZaO 



7ra; 

a .,-1 1 f°° z- 



x C(Va) - 1 ™~ + - / -r^dz , (7) 

where [z] is the integer part of z. Application of standard summation formulas ^| shows 
that, if a = 1, M^ 1} (0) grows logarithmically, whereas, if 1 < a < 2, 

M»(0) ~ (8) 

as n — > oo; finally in case of classical random walk (a > 2), M n (0) = (^(n 1 / 2 ) ^(|. Since 
flights are symmetric and start from 0, Pu is, for every k, an even function with a maximum 
in Q and hence M«(0) > M<» W , fa eve ry n and for eve ry , , 0; therefore, if a < !, 
.M^(x) < oo. A series expansion of Eq. (£Q) shows that 

M,l 1) (x)=M ? l 1 )(0)-l 

oo -p (2l±l\ 21 n 

♦Em^^E^ <») 

i=l v 7 k=l 

now for every < a < 2 and every I, 21 + 1/a > 1, and lim^oo X]fc=i A; s~ = C((2Z + l)/a) 
is finite. Then M ( n\x) = O(m£\0)), that is the last term on the RHS of © just takes 
into account the delay with which the flyer reaches site x; in particular, if a > 1, M n (0) 
diverges and Mn (x) ~ M n (0). In conclusion, a one- dimensional flight is transient if and 
only if a < 1, a result which has been obtained in a somehow more complex way in j^J. 

Consider now a (^-dimensional lattice, with d > 2, and assume the probabilities 
along the different coordinates to be independent; then Eq. (0 ) becomes Mn (0) = 
1 + (^T^ELi Note that > for < a < 2 and c/ > 2, the condition d/a > 1 

holds and hence M {d) {0) = 1 + (^^) ((d/a) is finite; M^ d) (0) as a function of n can be 
computed by using a method similar the one-dimensional case, and the result is that the 
trend is given by F(n) = ( j-^^- J yi 9 ^ — 1 j + 0(n~a). Finally it should be observed 
that the results for Mn\x), x ^ 0, obtained above, can be extended in a straightforward 
way to multidimensional lattices. Thus Levy flights on lattices of dimensions higher than 
one are always transient; if a > 2, Mn\o) = 0(log(n)) and, if d > 2, M n ' (0) converges to 
a finite value jl(3], and the walk is transient 
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Note that, when a — 1, M„ (0) has the same trend of M^(0) in the Gaussian regime, 
an instance of Levy flights increasing the effective dimension of the walk 

III. FINITE INTERVALS WITH ABSORBING BOUNDARIES 

In case of flights on a bounded set it is obvious that for reflecting boundaries 
diverges as n increases, since asymptotically pf; ~ 1/|£|> where \C\ is the number of sites 
whereas if boundaries are absorbing M.^ exists; here we shall consider just the case 
of one-dimensional lattices with absorbing boundaries. The map can be computed 

by means of numerical or analytical methods. In fact, Eq. (J2J) can be seen as a recursive 
method to compute Pf; and application of provides the result; alternatively, one can 
use the diffusion approximation to derive an analytical formula. Both methods have been 
used here and their results have been compared with S(x), the "experimental" number of 
visits generated by a Monte Carlo simulation. 

In a closed interval [—a, a] Eq. (j2J) becomes 

a 

p£ti(*)= X>(*-v)^ cl) (i/). ( 10 ) 

x=— a 

here, for reason of simplicity, instead of Q, we have used the transition probability, 
defined on integers y, 

p(y) = ^\y\- {a+1) if v^o, (n) 

and p(0) = 0, Z being a normalising constant. A similar form of p has been used in a 
work on the average time spent by flights in a closed interval In case of numerical 

calculations, obviously, the absolute length \y\ of a step must be truncated to some finite 
value: here max(|y|) = 2a, to allow flights to encompass the whole interval, and consequently 
Z = Y2y°=-2a M _( ' a+1 ' ) 5 V 7^ 0- Equation (fTTj) provides a valid transition probability for any 
a > and hence it can be used to model also classical Brownian motion; for a — ► oo 
the process becomes the simple symmetric walk. Note that by combining (J2J and (0) a 
recursive formula for Mn can be derived, namely M^^x) = YT x =-aPi x ~y)Mn (y) +<5 Xj o; 
however the separate use of and ((TJ) is to be preferred, in that it also yields values of the 
probability distribution and this is useful to check the correctness of the results. 

In the classical theory of random walk the diffusion approximation allows to replace 
Pjp(x) with the pdf P^(x,t), solution of the diffusion equation [l6|; analogously for Levy 



0, In 



flights a super diffusion equation can be derived (see, among others, |l2j |15|, [121), whose 
solution is a series of eigenf unctions fk of the operator V a [l5(. Setting P^'(x, 0) = 6(x — 0), 
the pdf is P^>(x,t) = / m (0)/ m (a;) exp(A m t). Define, in analogy with the discrete case, 

M (1) {x,t) = [ P®(x,T)dT, (12) 
Jo 

then M<V(x,t) = YZ,=\ ^m 1 /m(0)/ m (x)(exp(A m t) - 1) where \ k are the eigenval- 
ues of V a ; obviously, < 0, for all k, and the asymptotic formula is Ai^(x) = 
E^=ijA»l _1 /™(0)/»(x). 



In 



171 ] a solution P^(x,t) of the super diffusion equation has been presented that, for 



symmetric flights, is 



P {1) (x,t) = jJ2^P[-D a ^m/L) a t] 



L 

m=l 



x sm 



/ mii(x + a)\ . /mna\ 



here L = 2a is the length of the interval and D a the diffusion coefficient; application of Eq. 
(H2J) to (fT3j). with t — >■ oo, provides an explicit form for Ai^(x), 



Ml(x) = L^(^r^ sm {-^—) sm {—>- (14) 

m=l v ' x ' 

Calculations of A4^\x) from Eq. (jl4|) need the numerical value of the diffusion coefficient 
D a , and it can be derived from the average time T a flyer spends in the interval, related to 
D a by the formula 17] 

4 (-ir . n ,s 

ttA, \tt J (2m+l) a + v 1 J 

x ' m=l v ' 

since T is defined as T = J^ a dx J* °° P(x,t)dt = J^ a Ai ( - 1 \x)dx the approximation T « 
^" = _ a iS(x) can be used to obtain the numerical value of D a . 

Figures ((TJ and (j2J) show .M^ for a = 0.8 and a = 1.8 respectively. It can be seen that 
the graph of M. tends to a triangular shape as a increases; indeed for simple symmetric 



random walk, (a = oo, p(x) = l/26\ x \ t i) M.^ l \x) = a — \x\ 

Figure (JHJ) presents the graph of A4^(0) as a function of a; note that the inflection point of 
the curve occurs at a = 2, that is at the boundary between Levy flights and classical random 
walks. In other words, A4^(0) shows a "phase transition" from Levy flights, characterised 
by small number of visits, to the Gaussian regime where visits are more frequent. 



t ■ I 
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FIG. 1: Graphs of , in case of absorbing boundaries: points denote the Monte Carlo simulation, 
the dashed line the numerical method via eqs.^JandQ an d the full line the analytical solution 
(see !14j) . The parameters are a = 0.8, L = 51; The Monte Carlo simulation comprises 10000 trials, 
in Eq. jl|n = 2000, and the index m in Eq. (|14jl ranges from 1 to 20 
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FIG. 2: The same as Fig.fflbut a=1.8 

IV. CONCLUSION 

The results of this note clarify how the mean number of times a site is visited by a 
random flyer depends on the dimensionality of the lattice, the value of a and the boundary 
conditions. In particular, it has been shown that unrestricted Levy flights are always 
transient, but for the unidimensional case with a > 1; restricted flights are transient if the 
boundaries are absorbing. In the last case computations show that the direct numerical 
method agrees very closely with "experimental data" generated by the Monte Carlo 
simulation, whereas the agreement is worse for Eq. (|14)1. especially when a is small (see 
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FIG. 3: Graphs of MW(0), as a function of a, starting from a = 0.2, in case of absorbing 
boundaries. Points represent the Monte Carlo simulation, the dashed line the numerical solution 
of El and and the continuous line results from Eq. (|14j) . with a < 2. The horizontal line is the 
result for a = oo 

Figs. Q and 0J); this is not surprising, since Eq. fllOl) deals directly with discrete variables, 
whereas Eq. f)14j) results from the diffusion approximation. On the other hand, obviously, 
Eq. (jl4j) provides a more general, analytical formula for and not just a set numerical 
values. 

We thank the two anonymous referees for useful advice and criticism. 
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